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I. INTRODUCTION 

Einstein's equation determines geometries; hence its 
solutions are equivalent classes under space-time difFeo- 
morphisms of metric tensors. It is this invariance, how- 
ever, which imposes a particular aftermath on every ini- 
tial value formulation for Einstein's equation. The ge- 
ometrical equation must be first converted into a sys- 
tem having a well posed Cauchy problem, and so with- 
out the diffeomorphism invariance. A preferred foliation 
of spacelike hypersurfaces on the space-time is usually 
introduced in order that adapted coordinates break this 
invariance. Einstein's equation is then decomposed into 
constraint equations on the foliation hypersurfaces and 
evolution equations. While the constraints are uniquely 
determined by this procedure, the evolution equations 
are not. Some of these evolution equations turn out to 
be hyperbolic. This is in accordance with a main aspect 
of general relativity, that of causal propagation of the 
gravitational field. 

Hyperbolicity refers to algebraic conditions on the 
principal part of the equations which imply well posed- 
ness for the Cauchy problem, that is, the existence of 
a unique continuous map between solutions and initial 
data. There are several notions of hyperbolicity, which 
are related to different algebraic conditions. Some no- 
tions imply well posedness for the Cauchy problem in 
constant coefficient equations but not in more general 
systems, such as quasilinear equations. See 0, 01 for re- 
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views intended for researchers on general relativity. Re- 
garding quasilinear systems, strong hyperbolicity is one 
of the more general notions of hyperbolicity that implies 
well posedness for the Cauchy problem. The proof in- 
volves pseudodifferential analysis 0, ■ Symmetric hy- 
perbolic systems are a particular case of strongly hyper- 
bolic systems where well posedness can be proved without 
using pseudodifferential techniques. Several equations 
from physics can be cast into this symmetric hyperbolic 
form 5]. Finally, weak hyperbolicity is a less rigid no- 
tion than strong hyperbolicity but it does not imply well 
posedness for quasilinear equations. 

A definition of strong hyperbolicity for pseudodifferen- 
tial first order systems is introduced in Sec. |nj Differ- 
ential first order strongly hyperbolic systems known in 
the literature are included. A reason for this definition 
is that it incorporates precisely the hypothesis needed 
for proving well posedness. The proof involves standard 
pseudodifferential techniques. If an m-order differential 
system has a first order, differential or pseudodifferential, 
strongly hyperbolic reduction, then it is well posed. See 
the end of this section for an example of a first order pseu- 
dodifferential reduction of the wave equation. Also see 
Appendix ^ for a brief and self-contained introduction 
to the subject of pseudodifferential operators. 

Although strongly hyperbolic systems are at the core of 
the various proofs of well posedness for the Cauchy prob- 
lem in general relativity, they have played, until recently, 
no similar role in numerical relativity |6||. Finite differ- 
ence schemes have been implemented for non-strongly- 
hyperbolic equations. However, Lax's equivalence theo- 
rem does not hold in these situations i2|- It has been 
shown that discretization schemes standard in numeri- 
cal relativity are not convergent when applied to weakly 
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hyperbolic and ill posed systems H, • As more com- 
plicated situations are studied numerically, the interest 
in strongly hyperbolic reductions of Einstein's equation 
is increasing. There is also much experience and a vast 
literature in numerical schemes based on well posed for- 
mulations coming from inviscid hydrodynamics pH . Illj | . 
This experience can be transferred into numerical rela- 
tivity when strongly hyperbolic reductions are used. 

Early numerical schemes to solve Einstein's equation 
were based on variants of the Arnowitt-Deser-Misner 
(ADM) decomposition 12]. Only recently has it been 
proven that a first order differential reduction of the 
ADM evolution equations is weakly hyperbolic . This 
is the reason for some of the instabilities observed in 
ADM-based numerical schemes 0,13 ^3- In the first part 
of this work the ADM evolution equations are reviewed. 
A densitized lapse function is introduced, with the den- 
sity exponent held as a free parameter. The evolution 
equations are reduced to first order in a pseudodifferen- 
tial way. It is found that the resulting system is weakly 
hyperbolic for every prescription of a positive density 
lapse and a spacelike shift vector. This is summarized 
in Theorem ^ pseudodifferential techniques and mode 
decomposition are at the core of a proof directed to com- 
puting the eigenvalues and eigenvectors of the principal 
symbol of the evolution equations, and to check that the 
eigenvectors do not span the whole eigenspace. The mode 
decomposition also helps to understand why the addition 
of the Hamiltonian constraint into the system does not 
produce a strongly hyperbolic system. 

Baumgarte-Shapiro-Shibata-Nakamura-(BSSN-)type 
systems are introduced in the second part of this work. 
They are essentially the densitized ADM evolution equa- 
tions where some combination of connection coefficients 
of the three-metric is introduced as a new variable. The 
mode decomposition of the densitized ADM evolution 
equations and previous work on the linearized ADM 
equations suggest the introduction of this variable. 
It turns out to be related with the variable F* of 
the BSSN system, defined by Eq. (21) in 14]. (See 
also Il9l-) Similar variables have been introduced in 
[Tl[i8nia.l20ll2l|. Their evolution equation is obtained, 
as in the BSSN system, from commuting derivatives 
and then adding the momentum constraint. It is shown 
here that the addition of the momentum constraint 
transforms a weakly hyperbolic system into a strongly 
hyperbolic one. This is the main result of this work, 
and it is presented in Theorem |21 The first part of 
the proof follows the previous one for the densitized 
ADM evolution equations. Once the eigenvectors are 
computed, and it is verified that they do span the whole 
eigenspace, the proof continues with the construction of 
the symmetrizer. This construction is carried out with 
the eigenvectors. Finally, the smooth properties of the 
symmetrizer are verified. 

The hyperbolicity of a family of BSSN-type evolution 
equations has previously been studied with a different 
technique p3 ] . The equations were reduced to a differen- 



tial system of first order in time and in space deriva- 
tives. The lapse was densitized, the eigenvalues and 
eigenvectors of the principal symbol were computed, and 
it was verified that the latter do span their eigenspace. 
A smooth symmetrizer was computed for a subfamily of 
systems, showing strong hyperbolicity in this case. 

All notions of hyperbolicity mentioned here require 
rewriting the evolution equations as a first order system. 
This can be done in a differential or pseudodifferential 
way. Some pseudodifferential reductions to first order 
have the advantage that no extra equations are added 
into the system, so there are no extra constraints. This 
reduces the algebra needed to compute the symmetrizer. 
These techniques are well known in the field of pseu- 
do-differential calculus. They were first used in general 
relativity in ^15] , where linearized ADM evolution equa- 
tions were proved to be weakly hyperbolic. The example 
below presents the wave equation as a toy model to un- 
derstand how the pseudodifferential first order reduction 
works. See Appendix IXI and also Sec. 5.3 in Q, for 
other possible first order reductions. Consider the wave 
equation on for a function h, written as a first order 
system in time, in appropriate coordinates, that is, 

dth = k, dtk = A/i, 

with t e [0, cx)), e R^, and A = S^-^d^dj the flat Lapla- 
cian. Here (5'^ = diag(l,l,l). Fourier transform the 
system in x', 

dth^k, dtk^^\ufh, (1) 

where h{t,u}i) is the Fourier transform of h{t,x^) in the 
space variables, as defined in Appendix IA 21 The function 
k is defined in an analogous way, and ]w]^ = S^^ujiOjj. The 
key step is to rewrite Eqs. |^ as a first order system by 
introducing the unknown £ := i\uj\h, where i\uj\ is the 
symbol of the pseudodifferential operator square root of 
the Laplacian. One gets 

dt£ — i\uj\k^ dtk — i\uj\i. 

The system so obtained is a reduction to first order of 
the original first order in time wave equation. Notice 
that there is no increase in the number of unknowns, 
just a replacement of h by £, and correspondingly no ex- 
tra constraints are introduced. For the wave equation 
the result is a symmetric hyperbolic system. In the case 
of the ADM equations the resulting system is weakly hy- 
perbolic, with or without densitizing the lapse function, 
while for the BSSN-type equations one gets a strongly 
hyperbolic system. 

In Sec. |n] a precise definition for well posedness is 
introduced for first order quasilinear pseudodifferential 
systems. Strongly hyperbolic systems are also defined 
and the main theorem asserting well posedness for these 
systems is reviewed. Section IIII Al is dedicated to re- 
viewing the densitized ADM equations. The main result 
here is Theorem ^ asserting that the resulting evolution 
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equations are weakly hyperbolic. The role of adding the 
Hamiltonian constraint is briefly discussed. Section llll ij| 
is dedicated to introducing the BSSN-type system modi- 
fying the densitized ADM equations. The main result of 
this work, Theorem asserts that this BSSN-type sys- 
tem is strongly hyperbolic for some choices of the free 
parameters. The key point is the introduction of the mo- 
mentum constraint into the evolution equations. Section 
IIVI summarizes these results briefly. Appendix ^ is an 
introduction to pseudodiffcrcntial calculus. It summa- 
rizes the main ideas and highlights the main results. It is 
intended for physicists interested in learning the subject. 
It provides all the background knowledge to follow the 
calculations presented in this work. A summary of this 
type of pseudodifferential calculus was not found by the 
authors in the specialized literature. 

II. WELL POSEDNESS 

Hadamard first introduced the concept of well posed- 
ness for a Cauchy problem. It essentially says that a well 
posed problem should have a solution, that this solution 
should be unique, and that it should depend continuously 
on the data of the problem. The first two requirements 
are clear, but the last one needs additional specifications. 
First, there is no unique way to prescribe this notion of 
continuity. Although a topological space is all that is 
needed to introduce it, Banach spaces are present in most 
definitions of well posedness. This refinement simplifies 
the analysis while still including a large class of problems. 
Second, as nonlinear systems lack, in general, global in 
time solutions, one can at most expect a local in time 
notion of well posedness. Discussions on well posedness 
can be found in and See for a summary. 

This section is dedicated to reviewing the minimum set 
of definitions and results on well posedness for quasilin- 
ear pseudodifferential strongly hyperbolic systems which 
are needed to describe the equations coming from general 
relativity. It assumes that the reader is acquainted with 
the notions from functional analysis and pseudodifferen- 
tial calculus given in Appendix IXI 

Consider the Cauchy problem for a quasilinear first 
order pseudodifferential system 

dtu = p{t,x,u,dx)u, u\t=o^f, (2) 

where u, / are m-dimensional vector valued functions, 
TO > 1, and X represents Cartesian coordinates in M", 
n > 1. Here p{t,x,v,dx) is a smooth family of pseu- 
dodifferential operators in ip^p parametrized by i € M"*" 
and V G M™. Let p{t,x,v,Lu) and p^{t,x,v,Lu) be their 
symbols and principal symbols, respectively. 

If p is a differential operator with analytic coefficients, 
then the Cauchy-Kowalewski theorem asserts that there 
exists a unique solution for every analytic data /. How- 
ever, solutions corresponding to smooth data behave very 
differently depending on the type of operator p. For ex- 
ample, write the flat Laplace equation in R"+^ and the 



flat wave equation in K"+-'^ as first order systems in the 
form dtu = A^diU. The matrices are skew symmetric 
for the Laplacian, and symmetric for the wave operator. 
Therefore, solutions of the form u{t, x) = M(i)e'"'^ for 
the corresponding Cauchy problems behave very differ- 
ently at the high frequency limit. The solutions of the 
Cauchy problem for the Laplace equation diverge in the 
limit |w| — > oo, while the solutions of the wave equation 
do not diverge in that limit. (See footnote [13 ■) 

Let i3(E") be a Banach space with norm || ||, whose 
elements are vector valued functions from to MJ" . The 
Cauchy problem 0j is well posed in _B(R") if given initial 
data f{x) € i3(R") there exists a solution u(t,x) which 
is unique in i?(M") for each t € [0,T), for some T > 0; 
and given any number e > there exists 5 > Q such that, 
for every data f(x) G _B(R") satisfying ||/ — /|| < 5 there 
exists a unique solution u(t,x) G i?(M") for t G [0,T) x 
K" for some T > 0, with \T — T\ < e, and satisfying 
\\u{t) - u{t)\\ < e, for all t G [0, min(T', T)). This means 
that the solution depends continuously on the data in the 
norm || || . 

Well posedness is essentially a statement about the be- 
havior of the solutions of a Cauchy problem under high 
frequency perturbations of the initial data. Here is where 
pseudodifferential calculus is most useful to study solu- 
tions of the Cauchy problem. The high frequency part 
of the solution can be determined by studying the higher 
order terms in the asymptotic expansion of symbols. 

A wide class of operators with well posed Cauchy prob- 
lem is called strongly hyperbolic. A first order pseu- 
dodifferential system 0) is strongly hyperbolic if p G 
and the principal symbol is symmetrizable. This means 
that there exists a positive definite, Hermitian operator 
H(t,x,uj) homogeneous of degree zero in lu, smooth in all 
its arguments for uj ^ 0, such that 

where p* is the adjoint of the principal symbol p_^. 

The definition summarizes all the hypotheses on quasi- 
linear systems needed to prove well posedness. It is the 
definition given in Sec. 3.3.1 in |10|| for linear variable co- 
efficient systems, and the so-called symmetrizable quasi- 
linear systems given in Sec. 5.2 in 

Consider first order differential systems of the form 
dfU = A^{t, x)diU -(- B{t, x)u. The symbol is p{t, x, uj) = 
iA^{t,x)LUj + B{t,x), and the principal symbol is 
p_^{t,x,Lu) = iA^ujj. If the matrices A* are all symmet- 
ric, then the system is called symmetric hyperbolic. The 
symmetrizer H is the identity, and p^+p* — 0. The wave 
equation on a fixed background, written as a first order 
system is an example of a symmetric hyperbolic system. 
Well posedness for symmetric hyperbolic systems can be 
shown without pseudodifferential calculus. The basic en- 
ergy estimate can be obtained by integration by parts in 
space-time. 

If the matrices are symmetrizable, then the differ- 
ential system is called strongly hyperbolic. The sym- 
metrizer H — H(t, X, Lu) is assumed to depend smoothly 
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on uj. Every symmetric hyperbolic system is strongly 
hyperbolic. pseudodifFerential calculus must be used to 
show well posedness for variable coefficient strongly hy- 
perbolic systems that are not symmetric hyperbolic [Sj. 
The definition given two paragraphs above is more gen- 
eral because the symbol does not need to be a polyno- 
mial in UJ. The definition given above includes first order 
pseudodifferential reductions of second order differential 
systems. These type of reductions are performed with 
operators like A, A, or ^, defined in the Appendix 1X1 

In the particular case of constant coefficient systems 
there exists in the literature a more general definition of 
strong hyperbolicity 0, The principal symbol 

must have only imaginary eigenvalues, and a complete 
set of linearly independent eigenvectors. The latter must 
be uniformly linear independent in w 7^ over the whole 
integration region. Kreiss's matrix theorem (see Sec. 2.3 
in [lOj) says that this definition is equivalent to the ex- 
istence of a symmctrizcr H. Nothing is known about 
the smoothness of H with respect to i, x, and lo. The 
existence of this symmetrizer is equivalent to well posed- 
ness for constant coefficient systems. However, the proof 
of well posedness for variable coefficient and quasilinear 
systems does require the smoothness of the symmetrizer. 
There are examples showing that this smoothness does 
not follow from the previous hypothesis on eigenvalues 
and eigenvectors of . Because it is not known what ad- 
ditional hypothesis on the latter could imply this smooth- 
ness, one has to include it into the definition of strong 
hyperbolicity for nonconstant coefficient systems. 

A more fragile notion of hyperbolicity is called weak 
hyperbolicity, where the operator has imaginary eigen- 
values, but nothing is required of its eigenvectors. Quasi- 
linear weakly hyperbolic systems are not well posed. The 
following example gives an idea of the problem. The 2x2 
system dtu = Ad^u with a: G M and 




is weakly hyperbolic. Plane wave solutions of the form 
u{t,x) = u{t)e"^-'' satisfy \u{t)\ < \u{0)\{l + \uj\t). There- 
fore, plane wave solutions to a weakly hyperbolic system 
do not diverge exponentially in the high frequency limit 
(as in the case of Cauchy problem for the Laplace equa- 
tion) but only polynomically. This divergence causes so- 
lutions to variable coefficient weak hyperbolic systems to 
be unstable under perturbations in the lower order terms 
of the operator, as well as in the initial data. 

The main theorem about well posedness for strongly 
hyperbolic systems is the following. The Cauchy problem 
0) for a strongly hyperbolic system is well posed with 
respect to the Sobolev norm |j ||s with s > n/2 -\- 1. The 
solution belongs to C([0, T), H"^), and T > depends only 
on ll/lls- 

In the case of strongly hyperbolic differential systems, 
this is Theorem 5.2.D in The proof for pseudodiffer- 
ential strongly hyperbolic systems is essentially the same. 
One builds an estimate for the solution in a norm, defined 



using the symmetrizer, equivalent to the Sobolev norm 
H'^. Then the argument follows the standard proof for 
differential systems. The construction of the symmetrizer 
is basically the one carried out in ^2311. 

III. ADM DECOMPOSITION OF EINSTEIN'S 
EQUATION 

The ADM decomposition of Einstein's equation is re- 
viewed. The densitized lapse function is introduced in 
Sec. nil Al Theorem ^ says that the resulting evolution 
equations are weakly hyperbolic, for every choice of a 
positive densitized lapse function and spacelike shift vec- 
tor as given fields. The BSSN-type system is introduced 
in Sec. IIII Bl It is essentially the system given in 
without conformal-traceless decomposition, keeping the 
three connection functions and densitizing the lapse func- 
tion. It is reduced to a first order pseudodifferential sys- 
tem like the ADM evolution equations. The main re- 
sult, Theorem El asserts that BSSN-type equations are 
strongly hyperbolic, for every positive densitized lapse 
function and spacelike shift vector. 

Let {M,gab) be a space-time solution of Einstein's 
equation. That is a four-dimensional, smooth, orientable 
manifold M , and a smooth, Lorentzian metric gab solu- 
tion of 

Gab ~ f^Tabj 

with Gab = Rab — RQab/'^ the Einstein tensor. Tab the 
stress-energy tensor, and n = Stt. Ricci's tensor is Rab 
and R denotes Ricci's scalar. Latin indices a, 6, c, d 
denote abstract indices, and they are raised and lowered 
with g""^ and gabi respectively, with gacd'^^ = 

^ The 

unique torsion- free metric connection is denoted by V a- 
The conventions throughout this work are 2V[aV(,]Uc — 
Rabc^Vd for Riemann's tensor and (— , -h) for the 
metric signature. 

Prescribe on M a foliation of spacelike hypersurfaces 
by introducing a time function i, which is a scalar func- 
tion satisfying the condition that V at is everywhere time- 
like. Denote the foliation by St, and by Ua the unit 
normal to St such that n" — g°'^ni, is future directed. 
Therefore, Ua = —NWat for some positive function A'^. 
Fixing the foliation determines its first and second fun- 
damental forms hab — gab + nanb and kab = —ha'^'^cn-b, 
respectively. Decompose Einstein's equation into evolu- 
tion equations Q , Q , and constraint equations (jSl , lO , 
as follows 

^nhab = '2'kabi (3) 

^nkab — ^ ^ Rab 2/Ca kf^c -\- k kab 

-{DaDbN)/N~KSab, (4) 

Dbka'' - Dak = nja, (5) 
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(6) 



where denotes the Lie derivative along n", and k = 
ka°'. Here '^^^Rat, ^^^R, and Da are, respectively, the Ricci 
tensor, the Ricci scalar, and the Levi-Civita connection 
of hab, while h""^ denotes its inverse. The stress-energy 
tensor is decomposed as Sab = {ha'^hb'^Tcd — Thab/^) , 
with T = Tabg''\ Ja = -ha'n^'T.d, and p = Tabu'^nK 

Introduce on M a future-directed timelike vector field 
t°. Impose the additional condition f^Vat = 1, that is, 
the integral lines of are parametrized precisely by t. 
This condition implies that the orthogonal decomposi- 
tion of with respect to St has the form t°- — Nn°- -\- 
with na/?° = 0. is called the lapse function and /J" 
the shift vector. The integral lines of t"" determine a dif- 
feomorphism among the hypersurfaces St- This, in turn 
determines a coordinate system on M from a coordinate 
system on 5'o. Lie derivatives with respect to can be 
rewritten in terms of and fS"" . The resulting equations 
are called the ADM decomposition of Einstein's equation. 



A. Densitized ADM equations 

Consider the ADM decomposition of Einstein's equa- 
tion. Let be a coordinate system adapted to the folia- 
tion St, where — t and are intrinsic coordinates on 
each St that remain constant along the integral lines of 
i". Greek indices take values 0, 1, 2, 3, and latin indices i, 
j, fc, /, take values 1,2,3. In these coordinates, = Sq'^, 
= -7V(5^°; then = (5//3* and = {Sq^" - /3^)/A^. 
The components of the space-time metric have the form 

In these coordinates 



(3) 



Rij = ^h'''^ [-dkdihij - didjhki + 2dkd(^^hj)i] 

hi hi 

+jzknj - In iki , 



where ^^^'^ = /i^^ h^,^ hai^T^^ii," are the spatial com- 
ponents of the Christoffel symbols of g^^i,, and 7,;^*^ := 

Densitize the lapse function, that is, write it as = 
[HifQ, where h := det(/iij), & is a constant, and Q is 
a given, positive function. This modifies the principal 
part of Eq. |0J. New terms containing second spatial 
derivatives of hij come from {DiDjN)/N. 

Summarizing, the unknowns for the densitized ADM 
equations are hij and kij. The evolution equations, Eqs. 
©, igjl, have the form 



C(^t-p)hj = {N/2)h'''[-dkdiKj^{l + h)d^d,hk 
+2dkd(ihj)i] + B,j, 



(7) 



(8) 



where C(t-i3)hij = dtKj - {P'^dkhij + 2hk{idj)l3^), and 
the same holds for kij. The nonprincipal part terms are 
grouped in 

B^j = N [7,fc;7i''' - lij^lki - 2h'-kji 

-\-kijki Aij KS'yj , 
Aij = a,aj - j^j'^Qk + didj (In Q) + 2b'yiknj '^^^^ , 

with = n^V^n^ = Df^{\nN). The relations InA^ = 
b\nh+\nQ and {DiDjN)/N = {b/2)h''^d^djhki + Aij 
were used. 

The following result asserts that densitized ADM evo- 
lution equations are weakly hyperbolic. 

Theorem 1 Fix any positive function Q, a vector field 
/3' , and first and second fundamental forms hij , kij on 
Sq- lfb>0, then Eqs. Q), ^ are weakly hyperbolic. If 
b < 0, these equations are not hyperbolic. 

The proof has two steps: first, to write down Eqs. 10, 
as an appropriate first order pseudodifferential sys- 
tem, Eqs. ©, ifTUI) : second, to split the corresponding 
principal symbol into orthogonal parts with respect to 
the Fourier variable uji, and then to explicitly compute 
the associated eigenvalues and eigenvectors. 

Proof. First order reduction. Compute the symbol 
associated with the second order operator given by Eqs. 
0, ©, that is. 



dth,j= / {-2Nkj+iuJkl3''kj+2hi^,dj)(3'}e"^''duj, 

J St 



dthj = 



{{N/2) [\Lu\lh,, + (1 + b)LO,LO,h'''hkl 

't 

2uj^uj^ihj)k + iijJkP'^kij + Bij I e'^^'ticj. 



where hij and kij denote the Fourier transforms in .t* of 
hij and kij, and 



Bi] - B,j 



2ki^,d,)l3' 



denotes the terms not in the principal symbol. Here 
duj — dijj/(27r)'^/^, \ujW = uJiOjjK^^ , and we will use the 
convention lo"^ = ujjh^K Transform this second order 
symbol into a first order one via iij = i\Lo\fihij, where 
\uj\^^ = uJiUJjS^^ , with (5*^ — diag(l,l,l). The associated 
first order system is then 

Jtirj ^ [ {i\uj\h\-{2N/a)hj+pi,j 
JSt ^ 



+24(,9j)/5'}' 



(9) 



dth, = j [i\Lu\h [-iNa/2) (4 + (1 + b)Cd,Cojh'''iki 
-2u''Cji^Jj)k'j + (3kij + Bij^e^'^^'duj, (10) 
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with a = \u;\h/\uj\s, ui^ = LJi/\Lj\h, (3 := Wfc/?'", and % = 
i\uj\shije'^'^^duj. Then the symbol of equations above 
can be written as 

p{t,x,u,iuj) = i\u}\fiP-^{t,x,u,u}) + ]B{t,x,u,LLj), (11) 

where (IBu)'^ :— (2i'j.(i9j)/3'^, By ), ii^ :— {£ij, kij), with 
the upper index T meaning transpose. The principal 
part operator can be read out from the terms inside 
the square brackets in Eqs. ©, l|10(l . Notice that the 
definition of the principal symbol here differs from the 
one given in Sec. El by a factor of ilwl/i. (In particular, 
the eigenvalues of p^ as defined here must be real to be 
hyperbolic.) 

Eigenvalues and eigenvectors of p^ . Once the prin- 
cipal symbol is known, it only remains to compute its 
eigenvalues and eigenvectors. The assumption a = 1 fa- 
cilitates the computations. It is not a restriction since 
the norms | \s and | \h are equivalent and smoothly 
related, and therefore the properties of the eigenvalues 
and eigenvectors of the principal symbol are the same 
with either norm. Furthermore, one can check that if 
ii^ = {£ij, kij) is an eigenvector of ^^^(a = 1) with eigen- 
value A, then u^{a) = {a^^^'^lij , a^^'^ftij) is an eigenvec- 
tor of Pj^(a) with the same eigenvalue A. Therefore, from 
now on a = 1 is assumed. A second suggestion for doing 
these calculations is to decompose the eigenvalue equa- 
tion p_^u = Xu into orthogonal components with respect 
to uji. Introduce the splitting 

£ij = ujiujji + i'qij/2 + 2oj(ii'j^ + £'^^J^ , 



(12) 



k. 



ij = Lu^ibjk + k'qij/2 + 2w(ifc^.| -f- k'^^^^ , (13) 



where qij := hij — ujiujj is the orthogonal projector to (2;^, 
and 

£ = Cj'Cj'h^ , £' = g*J"4- , £[ = g>(I-'4i , 



The same definitions hold for the kij components. This 
decomposition implies that u = u^^^ + u'^' + u^^^ where 

,(1) ^ I" u}iLVj£+ [qij/2)£' 
uJiUjjk + {qij/2)k' 



(14) 





2u}(^i£'j') 




" f' 




. 2'I'(ifc'j) . 




V 



The principal symbol p_^ 



and the eigenvalue equation 
p^u = Xu can also be decomposed into the same three 
parts. The first part is four dimensional, corresponding 



to the variable u^^\ that is, the scalar fields, £, k, £', and 
fc'. The eigenvalues are 



±1, A 



(1) 



±Vb, 



where A := (A — I3)/N, so the role of the shift vector 
is to displace the value of the eigenvalue by an amount 
P — CokP'', and the lapse rescales it. But a change of lapse 
(which here is the function Q) and shift cannot change a 
real eigenvalue into an imaginary one. It cannot affect the 
hyperbolicity of the system. The associated eigenvectors 
for this first part are 



Al 



2[(1- 
T[(l 



(1 



b)q^J/2] 
-b)q^/2] 



2a;,: 
tV6 



uJiCjj 



Notice that for 6 = 1 the two eigenvectors m^^^^ collapse 

to the two eigenvectors u'^^ . The conclusion for this part 
is that the eigenvalues are real for 6 > 0, and the four 
eigenvectors are linearly independent for ^ 0, h ^ 1. 

The second part is also four dimensional and corre- 
sponds to the variable u'^'^\ that is, the vector fields £[ 
and k[. (The vector £[ has only two independent compo- 
nents because of the condition ^^w* = 0. The same holds 
for k'^.) The result is 



a(^) = 0, 



-(2) 
Al 



The eigenvalue has multiplicity 4, but there are only two 



linearly independent eigenvectors. Here, 



represent 



two linearly independent vectors, each one orthogonal to 
and labeled with the index which takes values 1, 2. 
This part is the main reason why the ADM equations are 
weakly hyperbolic. 

The last part is again four dimensional and corresponds 



to the variable u^^\ that is, the two-tensor fields £' 
and k 



1,^^^.^. (The tensor has only two independent 
components because of the symmetry, the orthogonality 
to and the trace- free condition. The same holds for 
k'l^ijy) The result is 



;(3) 



±1, 



^(3) 



The eigenvalues each have multiplicity 2, and there are 
four linearly independent eigenvectors. Here W(fc;)^ rep- 
resent two linearly independent symmetric, traceless ten- 
sors, orthogonal to a)^. 

At the end one gets the following picture. All eigen- 
values are real for 5 > 0. Notice that X^^ becomes imag- 
inary for 6 < 0, so the equations are not hyperbolic in 
this case. With respect to the eigenvectors, there are two 
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main cases. First, 6 > and b 1. Then, the eigenvec- 
tors of the first and third parts of do span their associ- 

ated eigenspaces; but the eigenvectors u\ corresponding 
to the second part of do not span their eigenspace. In 
the second case, 6 = or 6 = 1. In this case there are 
hnearly dependent eigenvectors even among the scalar 
variables. Therefore, the conclusion is that the system 

0, is weakly hyperbolic for 6 > 0. q 
It is interesting here to comment on the role of the 

Hamiltonian constraint. Suppose that a term of the form 
ahij times Eq. © is added to Eq. ^ . Here a is some real 
constant. Can this modification alter the hyperbolicity 
of the ADM equations? One might think that adding the 
Hamiltonian constraint to the ADM evolution equation 
could have a similar role as densitizing the lapse function 

1. e., it could keep both eigenvectors u^-^^ linearly indepen- 
dent. The fact is, it does not. Such an addition of the 
Hamiltonian constraint modifies only A^^^ and u^xj, and 
does not modify u^^^-* and x'^\ The result is 

„(i) _ r 2[{1 + b + 2a)u},Cjj + (1 ~b + 2a)qij/2] ' 
~ [ T\/TT2^[(1 -I- b)i:;jiQj + (l - b)q,j/2] _ • 

Therefore, adding the Hamiltonian constraint only helps 
to keep the eigenvectors u'"^^ independent of the so 
it helps only in the case 6=1, where the former collapse 
onto the latter (for a = 0). For b ^ 1 the addition of the 
Hamiltonian constraint does not contribute to make the 
vectors u^^^ linearly independent, whereas densitizing the 
lapse does. 

B. BSSN-type equations 

Consider the densitized ADM evolution equations lO, 
©. Introduce into these equations the new variable 

By definition n^/'^ = 0, that is, /" ~ 0, so the new vari- 
ables are the components /' = h^^ [h^'- dkhij — dj{lnh)], 
where h — yj det(/iy) as above. They are related to 
the three connection variables F* of the BSSN system 
defined in Eq. (21) in 0|. More precisely, '^ijt^ — 
fi + {l/3)di{\nii), where jij is defined in Eq. (10) of 
that reference. The evolution equation for fi is obtained 
by taking the trace in indices v, a of the identity 

+Kaa^ + —hf^sh,," he" du'da'P^, 

where CnluJ = n''d^j^J+2j^^Jd^-)n''-j^a''d^n^ , and 
adding to the result c times the momentum constraint 
(O. Here c is any real constant. One then gets 

CnU = (c - 2)i?,fc^'' + (1 - c)D^k + C^, 



where the nonprincipal terms are grouped in 

Ci_i = cKjfi 2k^i,a ~\- ka^j^ 2^i,„^k 
-2k^,r + {l/N)h^,,h'''d„ds(3r 

Summarizing, the unknowns for BSSN-type systems 
are hij , fcy , and fi . The evolution equations are 

^(t-/3)/iy = -27Vfcy, (15) 

'0-(t~f3)hj = ^h^^[-dkdihij -bdidjhki] 

+Nd^J^)+B,j, (16) 

Cit-p)f. = mc-2)h'''dkhj + {l~c)h''^d,kkj]+C,, (17) 

where C(t-|3)f^ = dtft-{P^ d-j fi+fjdif3^), while C(t-i3)htj 
and Ci^t-(3)kij are defined below Eqs.Q-lIHJl, and 

Bij = N [2^ki{il])^^ + likilj^^ - lijilk"^ 
2ki kji + kijki A.ij K/Sjjj , 
C, = A^[C, + (c-2)(7fe/fc,^-7fe,^fc/)]. 

The constraint equations are Eqs. Q, ^ and 

- /^""7.^" - 0. 

The main result of this work asserts that BSSN-type 
evolution equations are strongly hyperbolic for some 
choices of the free parameters. 

Theorem 2 Fix any positive function Q , vector field (3^ , 
first and second fundamental forms hij, kij on Sq. 

If b > 0, b ^ 1, and c> 0, then Eqs. C^-^I^ are 
strongly hyperbolic. 

Assume that b ^ I. If c ^ 2, then Eqs. < j5|) - |7^ are 
strongly hyperbolic; if c 2, c > 0, then they are weakly 
hyperbolic. 

One can check that the system (|15|) - (|17|l remains 
strongly hyperbolic under a transformation of the form 
Fi = fi + d di (In h) for any real constant d, in particu- 
lar d — 1/3, which gives the BSSN variable F*. (See the 
comment at the end of this section.) 

The first part of the proof follows the argument that es- 
tablishes Theorem^ That is, from Eqs. H15|) - (|17|l obtain 
the first order pseudodifferential system Eqs. H18|I - H2U|) 
below. Then compute the eigenvector and eigenvalues, by 
splitting the principal symbol into orthogonal parts with 
respect to the Fourier variable oji. Finally, the second 
part of the proof is the construction of the symmetrizer. 

Proof. First order reduction. Compute the symbol 
associated with the second order operator given by Eqs. 
(Cni-lIIZI): 

dth,j = j [-2Nkj + luJkP'^k, + 2/ifc(,a,)/3^-} e*""(?w. 
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dthj = 



{(iV/2) [\uj\lh,,+bu;,ujjh'''hki 



dtf. 



St 



|i7V (c - 2)kikUJ^ + (1 - cjuj^h'^^kkj 



where the terms not in the principal symbol have the 
forms 



Here du; = duj/{2'K)^/^, = LJiLjjh'^^ , and — ujjh^^ . 
Introduce the unknown £ij — i\uj\shij, with \uj\g = 
ujiUjjS'^^ , where S^^ — diag(l,l,l). The resulting pseu- 
dodifferential system is a first order one, given by 



(18) 



{i\uj\h [iNa/2) (-4 - b^.Cojh'^'iki 
+2Cj(^ J,)^ + Pho] + &y } e'^'^duj, (19) 



\i\u\h \N({c-2)kk^ 
+ (1 - c)u,h''^kkj) + (3^] + C,} e'^^du;(20) 



with a = \u!\h/\c^\s, i^i = ^^j/l^^l, ^ = and % = 

i|tj|5/iije*'^"^atj. The symbol of Eqs. (HHJl-lEnil has the 



form, 



p{t,x,u,iuj) = i\uj\hP-^^{t,x,u,Lu) + lB{t,x,u,Lu), (21) 



where (Bu)^ := (24(^9J)/?^ S„ , C,), 



) : f^ij 1 fi ): 



the index T denotes the transpose, and the principal sym- 
bol can be read out from the terms inside the square 
brackets in Eqs. (|18|I - H2UI) . As in the proof of Theorem 
^ the definition of the principal symbol here differs from 
the one given in Sec. El by a factor of ilujli^. 

Eigenvalues and eigenvectors of p^. Following the 
proof of Theorem ^ a = 1 is assumed. One can check 
in this case that, if = {(ij,kij,fi) is an eigenvec- 
tor of p^{a = 1) with eigenvalue A, then u^{a) — 

{a~'^^^£ij, a^^^kij, a^^'^ fi) is an eigenvector of p^(a) with 
the same eigenvalue A. 



The orthogonal decomposition (|12|l , H13|l simplifies the 
calculation. In addition, decompose 

h=Co,f + fi 

where / = ujif^ and /■ = qi^ fj. This decomposition 
implies that u — u^^^ 



vS^^ + iS^") where 



Cj,Cojl+{q,j/2)e 
ujiUjjk + {qij/2)k' 
i^if 



2Cj(J' 



2Cj(ik'j) 







(22) 



Split the principal symbol p^ and the eigenvalue equation 
p_^u — Xu into the same three parts. The first part is five 

dimensional, corresponding to the variable u^^\ that is, 
the scalar fields /, £, k, and k'. The eigenvalues are 



A«=±l, 



aW = ± A 



A 



(1) 



0, 



each having multiplicity 1, where again A := (A — $)/N. 
Hence, the eigenvalues are real if & > 0. The correspond- 
ing eigenvectors are 



2[{h^c+l)Cj,Uj + {\-h){q,jl2)] 
T[{h-c+l)i:^,Coj + (1 - h){q,j/2)] 
(2 — c)hCji 



2ujiUjj 




2ujiU)j 











_ (1 + b)Cj, _ 



Notice that both eigenvectors u^^^ collapse if 6 = 0. Also 

see that the eigenvectors u^^^ collapse to u^-^^ in the case 
b ~ 1 and 2. Thus, in these cases Eqs. ((T^ - (pnil are 
weakly hyperbolic. In the case 6 = 1 and c = 2 the eigen- 
values ±Vb collapse to ±1. Therefore, one hasA^^^ = ±1, 
each with multiplicity 2, and A2^'' — 0, with multiplicity 
1 . There are five linearly independent eigenvectors in this 
case, 



2ujiOjj 








Tqio/2 













The second part is six dimensional and corresponds to 
the variables u'-^\ that is, the vector fields i'^, and k'^, 
orthogonal to Coi. The eigenvalues are 



A 



(2) 



0, a(2)=±Vc/2, 
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where each one has multiplicity 2. They are real if c > 0. 
There are six linearly independent eigenvectors in the 
case c > 0, given by 



induces the same splitting in = p^'^"' +P^ 



and therefore in H = H^^^ + + ^(3)^ 
The result is 



rr (3) 



-(2) 






"(2) 
' "A. = 







(2 - c)v/ 



H 



TT (i)kl 



TT (i)kl 



W {f)k 



where Vj"^ represent two linearly independent vectors, 
each one orthogonal to (D^, labeled by the index A which 
takes values 1,2. Here is the key role of the momentum 
constraint. If c = 0, that is, the momentum constraint is 
not added to the system, then the two eigenvectors 
become linearly dependent, as occurs in the densitized 
ADM evolution equations. 

The last part is four dimensional and is the same as in 
Theorem n It corresponds to the variables u'^'^\ that is, 
the tensor fields ^'^ -^^ and fc'^ .^.^ . The result is 



where 



TT (i)kl 

TT {k)kl 

TT (e)ki 

TJ U)k 



ttW {i)kl 



tt{-2) (e)kl , r7-(3) {e)kl 



rr(l) {k)kl 
(k)ij 



r7-(2) (fe)fe; 
^{k)^j 



rr(l) {i)kl , tt{-2) {i)kl 

(in (/)« 



r7-(3) (fe)fe; 
^{k)^j 



x'? = ±1, 



.(3) 

= 







The eigenvalues each have multiplicity 2, and there 
are four linearly independent eigenvectors. The ten- 
sors represent two linearly independent symmetric, 
traceless tensors, orthogonals to uji. 

The symmetrizer. The operator H = (T^^)*(T^^) is a 
symmetrizer of system (|21f) . where T is an operator whose 
columns correspond to the eigenvectors of p^, is its 
inverse, and * denotes the adjoint. Then it only remains 
to do a lengthy, although straightforward, calculation. 

There are two hints that help to simplify the con- 
struction of the symmetrizer. They are based on the 
observation that the principal symbol has the particular 
form = Np_^ + [31 ^ where I is the identity matrix 

J = Aia,g{h^i'hj)\ hj^j^h^) ,hi'), and 



~ {i)kl Q 



(k)kl Q 



t^l{k)i3 

p^^^^ik)M 



if)k 



where the indices {£), (fc), and (/) indicate rows and 
columns, that is, equations and variables, respectively, 
and the matrix components are given by 



(k)kl _ 



-l{k)ij 



if)k 



t^iik)tj 



kl 



Then, the first hint is that a symmetrizer for p^ is indeed 
a symmetrizer for the nondiagonal elements p_^ . A second 
hint is that the orthogonal decomposition in Eq. I|22|) 



The matrix coefficients of each part depend on a);. The 
scalar variable part of the symmetrizer has the form 



o-(l) (l)kl 



H-(l) {f)k 

rr(l) {k)kl 
^(k)io 



ttW (l)kl 



lij „kl 



HkkCo^ujjU^u)'- + Hk'k' y-<7 
2 



where the scalar functions that appear above are the fol- 
lowing: 



Hu 



H 



ff 



Hfipi — 



_3_ 

262' 



Hk'k' — 



w 



462 
1 

462 



(6+l-c)2 + 6 
26(1 - 6)2 ' 

(6^ + 1 - cf + \? 
2(1-6)2 



(62 + l-c)(l + 6) 



Hfi' = — — 



Hkk' 



2(c-l) 



2(1-6) 
(62 + 1 - 



462 

jc-l-b) 
26(1-6) ' 



H 



(1-6) 



1 



-ip(^ + 3). 
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The symmetrizer for the vector variable part is 
2(c-2f + l 



f7-(2) {£)kl 

Tri'i) {k)kl 

^(2) (l)kl 
(f)i 

(/)« 



4c2 



4(c-2)-l , , 

— — ^m) . 

^2(i(,g,)('=(i'), 

4(c-2)-l n 



2c2 



9 

2^ 



Finally, the two-tensor part of the symmetrizer is 

r(3) (t)kl _ 1 



jj{3) (k)kl ^ 
(k)ii 



k I 



k I 
9(» ■ 



One can check that the symmetrizer so defined satisfies 
H{i\uj\iip^) + {i\uj\hP-^^)* H = and so trivially belongs to 
Si Q. This symmetrizer H{t, x, u, Co) is bounded for b ^ 1, 
c > 0, and smooth in all its arguments. The case 6=1 
and c = 2 can be computed in the same way described 
above, and the same conclusion holds. Then, for these 
two cases the system (|18|l - (|2Uf) is strongly hyperbolic, q 

Two further generalizations are immediate. The first 
one involves the Hamiltonian constraint. Suppose the 
term ahij times the Hamiltonian constraint @ is added 
to Eq. (|16|l . Here a is any real constant. What are 
the eigenvalues and eigenvectors of the resulting princi- 
pal symbol? The result is, as one expects, that the only 
change is in the scalar variable part p^^^^^. The eigenval- 



A^^^ =±y/l + 2a{2-c), A^- 



(1) 



each one having multiplicity 1, and A^^^ is real provided 
a(2 — c) > —(1/2). The corresponding eigenvectors are 



2[(6-c+l)(::;,tij- + (l-6)(7,y/2] 
TA[(& ~c+ l)Cj,Cbj + (1 - b)q.,j/2] 
(2 — cjbuji 



+a{2-c) 



2((D,;a)j + Qij) 
=FA(cJj(I'j + Qij) 
(2c- 1)lD, 



A3 



2uJiCjj 




2 [(3 — 2b)LJiUJj + bqij] 





+ a 





_ (1 + 6)cl;., _ 




3(Di 



The eigenvector uj^^ 
not change. The conclusion is summarized below. 



where A = ^Jl + 2a(2 - c). 



does 



Corollary 1 Consider Eqs. f J5)] - jJ7| ) and assume the 
hypothesis of Theorem\^ Assume b > and c > 0. As- 
sume that a term ahij times the Hamiltonian constraint 
^ is added to Eq. ^16]) . where a is a real constant sat- 
isfying a(2 — c) > —1/2. Then the resulting principal 
symbol, as defined in this Section, has real eigenvalues 
and a complete set of linearly independent eigenvectors. 

The second generalization involves the transformation 
Fi — fi-\-ddi{lR h), where d is any real constant. That is, 
instead of defining the BSSN-type system with the vari- 
able /i, define it with Ei. The new evolution equations 
have a different principal symbol from Eqs. I|15|) - (|17|) of 
the BSSN-type equations. However, one can check that 
it does not change the hyperbolicity of the system. In- 
deed, it modifies only the p^^*-^-* part of the principal sym- 
bol. Its eigenvalues remain the same, namely, A^^^* — ±1, 

Aa^"* — ±Vb, and Ag^"* = 0. The associated eigenvectors 
are now given by 



41) = 

Ai 



A2 



2[{b ^ c + 1)Cj,Cjj + (1 - b){q,j/2)] 
T[(6 - c + l)Cu^iJjj + (1 - fe)(g»j/2)] 
{2 - c){b + d)uj. 



2u)iLbj 




^^/buJiLJj 




_ (1 -1- d)u, _ 



2LijiLUj 





(1 + 6 + d)oJi 



Therefore, the hyperbolicity of the BSSN-type equations 
is not changed by this transformation. 



IV. DISCUSSION 

The first order pseudodifferential reduction performed 
in the space derivatives is here the main tool used to 
study the hyperbolicity of the BSSN-type systems. This 
technique is widely used in pseudodifferential analysis. It 
does not increase the number of equations, so there are 
no new constraints added to the system. It emphasizes 
that well posedness essentially captures the absence of 
divergent behavior in the high frequency limit of the so- 
lutions of a given system. This tool is applied to Eqs. 
I)15|l - H17|l . which have derivatives of first order in time 
and second order in space. They are obtained from the 
ADM equations by densitizing the lapse function and in- 
troducing the three connection variables fi. Its evolution 
equation is obtained by adding the momentum constraint 
to an identity from commuting derivatives. The positive- 
density lapse function and the spacelike shift vector are 
arbitrary given functions. There are free parameters 
given by the exponent in the densitized lapse and the 
factor in the addition of the momentum constraint. The 
resulting first order pseudodifferential system is strongly 
hyperbolic for some values of the free parameters. (See 
Theorem 0) 

The introduction of fi as a new variable is inspired by 
the variable of the BSSN system, defined by Eq. (21) 
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in 01 , and in the study of the linearized ADM evolution 
equations given in |15|. This variable is the crucial step 
that allows us to introduce the momentum constraint 
into the system. These two things, in turn, produce the 

(2) 

result that the vector variable eigenvectors u\ do span 
their eigenspace. This is the key feature that converts the 
weakly hyperbolic densitized ADM system into a strongly 
hyperbolic one. This property does not change when a 
term of the form d di (In h) is added to the system. Then, 
both results suggest why the BSSN system is preferred 
to the ADM equations for numerical analysis. This con- 
clusion agrees with a previous result in [22j. where the 
hyperbolicity of the BSSN system was also studied. It is 
shown there that a differential reduction to first order in 
time and space derivatives, together with a densitization 
of the lapse, produce a strongly hyperbolic system. The 
results in the present work are also consistent with nu- 
merical studies on the evolution equations presented in 
[l7. , 18 j 19 J . For further developments on these systems, 
see poll2l( . 

Finally, the role of the Hamiltonian constraint, when 
added to the ADM and BSSN- type evolution equations, 
is studied. In the case where the density lapse expo- 
nent 6 7^ 1 it does not affect the hyperbolicity properties 
of the two systems. That is, densitized ADM evolution 
equations remain weakly hyperbolic, and the BSSN-type 
system remain strongly hyperbolic. In the case 6=1 
the addition of the Hamiltonian constraint in both sys- 
tems prevents the two eigenvectors of the scalar variable 
block from collapsing onto each other. This keeps the 
BSSN-type equation strongly hyperbolic even in the case 
6 = 1, but is not enough to change the weakly hyperbolic 
character of the densitized ADM evolution equations. 
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APPENDIX A: ESSENTIALS OF 
PSEUDODIFFERENTIAL OPERATORS 

1. Introduction 

Pseudodifferential operators are a generalization of dif- 
ferential operators that make use of Fourier theory. The 
idea is to think of a differential operator acting upon a 
function as the inverse Fourier transform of a polynomial 
in the Fourier variable times the Fourier transform of the 
function. This integral representation leads to a gener- 
alization of differential operators, which correspond to 
functions other than polynomials in the Fourier variable, 
as long as the integral converges. 



In other words, given a smooth complex valued func- 
tion p{x, oj) from M" x M" with some asymptotic behavior 
at infinity, associate with it an operator p(a;, dx) : S ^ S. 
Here S is the Schwartz space, that is, the set of complex 
valued smooth functions in M", such that the function 
and every derivative decay faster than any polynomial at 
infinity. The association p ^ p, that is, functions into 
differential operators, is not unique. This is known to 
anyone acquainted with quantum mechanics. Different 
maps from functions p{x,ll!) into operators p{x,dx) give 
rise to different theories of pseudodifferential calculus. 
Every generalization must coincide in the following: The 
polynomial p{x,ijj) = X]|Q|<m where a is a 
multi-index in R" must be associated with the operator 
p{x,dx) = ^\a\<m'^oi{x)dx , that is, with a differential 
operator of order m. The map p ^ p used in these notes 
is introduced in subsection lA 31 of the Appendix. It is 
the most used definition of pseudodifferential operators 
in the literature, and the one most studied. 

The Fourier transform is used to rewrite the differ- 
ential operator because it maps derivatives into multi- 
plication, that is, [dxu{x)\'' = iLou{uj). This property is 
used to solve constant coefficient partial differential equa- 
tions (PDEs) by transforming the whole equation into 
an algebraic equation. This technique is not useful on 
variable coefficient PDEs, because of the inverse prop- 
erty, that is, = id,^u{uj). For example, one has 
[dxU + xu]" — i[d^u + ujii], and nothing has been sim- 
plified by the Fourier transform. That is why one looks 
for other ways of rewriting differential operators. The 
generalization to pseudodifferential operators is an ad- 
ditional consequence. Other transforms can be used to 
define different generalizations of differential operators. 
For example, Mellin transforms are used in p^. 

The functions p{x,oj) are called symbols. Differential 
operators correspond to polynomial symbols in uj. They 
contain the main equations from physics. Even strongly 
hyperbolic PDEs have polynomial symbols. Why should 
one consider more general symbols? Because the gen- 
eralization is evident, and it has proved worth doing it. 
The Atiyah and Singer index theorem is proved using 
pseudodifferential operators with smooth symbols, which 
are more suitable for studying homotopy invariants than 
polynomial symbols (25>] . Techniques to prove the well 
posedness of the Cauchy problem for a strongly hyper- 
bolic system require one to mollify polynomial symbols 
into smooth nonpolynomial ones The main applica- 
tion in these notes is simple: to reduce a second order 
partial differential equation to a first order system with- 
out adding new characteristics into the system. This is 
done by introducing the operator, a square root of the 
Laplacian, which is a first order pseudodifferential, but 
not differential, operator. The main idea for this type of 
reduction was introduced in [S^ . 

How far should this generalization be carried? In other 
words, how is the set of symbols that define the pseudod- 
ifferential operators determined? The answer depends 
on which properties of differential operators one wants 
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to be preserved by the general operators, and which ad- 
ditional properties one wants the latter to have. There 
are different spaces of symbols defined in the literature. 
Essentially all of them agree that the associated space 
of pseudodifferential operators is closed under taking the 
inverse. The inverse of a pseudodifferential operator is 
another pseudodifferential operator. This statement is 
not true for differential operators. The algebra devel- 
oped in studying pseudodifferential operators is useful to 
compute their inverses. This is important from a physical 
point of view, because the behavior of solutions of PDEs 
can be inferred from the inverse operator. One could even 
say that pseudodifferential operators were created in the 
middle 1950s from the procedure to compute paramet- 
rices to elliptic equations. [A parametrix is a function 
that differs from a solution of the equation p{u) = 5 by a 
smooth function, where p is a differential operator, and 
5 is Dirac's delta distribution.] To know the parametrix 
is essentially the same as to have the inverse operator. 

Spaces of pseudodifferential operators are usually de- 
fined to be closed under composition and transpose, and 
to act on distribution spaces and on Sobolev spaces. 
They can be invariant under diffeomorphism, so they can 
be defined on a manifold. This definition for pseudodif- 
ferential operators is not so simple as for differential op- 
erators, because the latter are local operators and the 
former are not. Pseudodifferential operators are pseu- 
dolocal. An operator, p, acting on a distribution, u, is 
local if p{u) is smooth in the same set where u is smooth. 
Pseudolocal means that the set where p{u) is smooth in- 
cludes the set where u is smooth. This means that p 
could smooth out a nonsmoothness of u. MoUifiers are 
an example of this kind of smoothing operator. They are 
integral operators, which justifies the name of pseudolo- 
cal. Differential operators with smooth coefficients are 
an example of local operators. The proofs of all these 
properties of pseudodifferential operators are essentially 
algebraic calculations on the symbols. One could say 
that the main practical advantage of pseudodifferential 
calculus is, precisely, turning differential problems into 
algebraic ones. 



2. Function spaces 

The study of the existence and uniqueness of solutions 
to PDEs, as well as the qualitative behavior of these so- 
lutions, is at the core of mathematical physics. Function 
spaces are the basic ground for carrying on this study. 
The mathematical structure needed is that of the Hilbert 
space, or Banach space, or Prechet space, which are com- 
plete vector spaces having, respectively, an inner prod- 
uct, a norm, and a particular metric constructed with a 
family of seminorms. Every Hilbert space is a Banach 
space, and every Banach space is a Frechet space. The 
main examples of Hilbert spaces are the space of square 
integrable fimctions L^, and the Sobolev spaces H'^, with 
k a positive integer, which consist of functions whose k 



derivatives belong to L^. The Fourier transform makes 
it possible to extend Sobolev spaces to real indices. This 
generalization in the idea of the derivative is essentially 
the same as one uses to construct pseudodifferential op- 
erators. Examples of Banach spaces are L^, spaces of 
p-power integrable functions, where L"^ is the particu- 
lar case p = 2. The main examples of Frechet spaces 
are C°°(f2), the set of smooth functions in any open set 
O C M", with a particular metric on it (the case il = R" 
is denoted C°°), the Schwartz space of smooth functions 
of rapid decrease, and its dual as a Prechet space, which 
is a space of distributions. 

This section presents only Sobolev spaces, first with 
non-negative integer index, and the generalization to a 
real index. The Fourier transform is needed to generalize 
the Sobolev spaces. Therefore, Schwartz spaces are intro- 
duced to facilitate the definition of the Fourier transform, 
and to extend it to L^. The next section is dedicated to 
introducing pseudodifferential operators. 

Let be the vector space of complex valued, square 
integrable functions on R ", that is, functions such that 
||u|| < oo, where \\u\\ := ^/ {u, u) and 

{u,v) := / u{x)v{x)dx, 

with u the complex conjugate of u. This set is a Hilbert 

space, that is, a complete vector space with inner prod- 
uct, where the inner product is given by ( , ) and is 
complete with respect to the associated norm || ||. 

The Sobolev spaces H'', for k a non- negative integer, 
are the elements of such that 

\\u\\l:= l|5"«f <oo, 

|a|<fc 

where a = (ai, a2, . . . , a„) is a multi- index, and for 
every such multi-index d°' denotes d"^ . . . d"'^ and 
|a| = Qij. The inner product in defines an inner 

product in iJ*^ given by 

|a|<fe 

Let S be the space of functions of rapid decrease, also 
called the Schwartz space, that is, the set of complex 
valued, smooth functions on K", satisfying 

\u\k,a ■■= sup 1(1 + |a;|^)''/^a"u| < cxd 

for every multi-index a, and all fc € N natural, with \x\ 
the Euclidean length in R". The Schwartz space is useful 
in several contexts. It is the appropriate space to intro- 
duce the Fourier transform. It is simple to check that 
the Fourier transform is well defined on elements in that 
space, in other words, the integral converges. It is also 
simple to check the main properties of the transformed 
function. More important is that the Fourier transform is 
an isomorphism between Schwartz spaces. As mentioned 



13 



earlier, the Schwartz space provided with an appropriate 
metric is an example of a Frechet space. Its dual space 
is the set of distributions, which generalizes the usual 
concept of functions. 

The Fourier transform of any function u e 5 is given 

by 

T[u][x) = u{x) := j e-"-'^u(w)(?tj, 

where 3lo = daj/(27r)"/^, while du and x-lu — SijX^u^ are 
the Euclidean volume element and scalar product in R", 
respectively. The map J- : S ^ S is an isomorphism. 
The inverse map is given by 

T~'^[u]{x) = u{x):= I e^'-'-'^u{uo)duj. 

An important property of the Fourier transform useful 
in PDF theory is the following: [d"u{x)]'' = i'°'^lj°'u{lo), 
and [x°'u{x)]'' = i^""^ d^u{uj) ^ that is, it converts smooth- 
ness of the function into decay properties of the trans- 
formed function, and vice versa. The Fourier transform 
is extended to an isomorphism J- : ^ , first proving 
Parseval's theorem, that is, (u, v) = (u, v) for all u, v ^ S 
(which gives Plancherel's formula for norms, ||m|| = 
in the case that the norm comes from an inner product, 
with u — v) and second recalling that S is dense in L^. 

The definition of Sobolev spaces for s real is based 
on Parseval's theorem. First recall that every u £ H*' 
with non-negative integer k satisfies d°'u € for |a| < fc, 
so Parseval's theorem implies \u)\^u{lo) £ L^. Second, 
notice that there exists a positive constant c such that 
{1/c){lu) < (1 + |w|) < c(w), where (lu) = (1 + |c^|2)i/2. 
Therefore, one arrives at the following definition. The 
Sobolev space i?^ for any s e M consists of locally square 
integrable functions in M" such that {bj)^u S LP'. This 
space is a Hilbert space with the inner product 

(u,w)s:= / {uj)'^'^u[ijj)v[ijj)duj^ 

and the associated norm is denoted by 

\\u\\l:= I M^^I^Hpdc.. 

One can check that H" C H'' whenever s' < s. Notice 
that negative indices are allowed. The elements of those 
spaces are distributions. Furthermore, the Hilbert space 
i/"" is the dual of i/*. Finally, two more spaces are 
needed later on, H-°° := UsmH'' and HseRH'. 
These spaces are, with appropriate metrics on them, 
Frechet spaces. A closer picture of the kind of element 
these spaces may contain is given by the following ob- 
servations. The Sobolev embedding lemma implies that 
H°° C C°°, while the opposite inclusion is not true. Also 
notice that S C H°°, and therefore H~°° C S' , so the 
elements of H^°^ are tempered distributions. 



3. Pseudodifferential operators 

Let S"™, with m G M, be the set of complex valued 
smooth functions p{x,uj) from R" x R", such that 

\d^dyix,co)\<c^{cor-\"\, (Ai) 

with Ca a constant depending on the multi-index a, and 
(lu) = (1 + Icjp)^/^. This is the space of functions whose 
elements are associated with operators. It is called the 
space of symbols, and its elements p{x, uj) symbols. There 
is no asymptotic behavior needed in the x variable, be- 
cause Fourier integrals are thought to be carried out in 
the u! variable. The asymptotic behavior of this variable 
is related to the order of the associated differential oper- 
ator, as one can shortly see in the definition of the map 
that associates functions p{x, uj) with operators dx)- 

One can check that S™ C S"™ whenever m' < m. Two 
more spaces are needed later on, 5*°° :— UmmS'^ and 
S °° := HrngMS"". 

Given any p{x,u!) G S"^, the associated operator 
p{x,dx) : S ^ S is said to belong to tp"^, and is de- 
termined by 

p{x,dx){u)^ ( e"->(x,w)u(w)cra;, (A2) 

for all u £ S. The constant m is called the order of the 
operator. It is clear that u £ S implies p{x, dx){u) £ C°°; 
however, the proof that p{u) e 5 is more involved. One 
has to show that p{u) and its derivatives decay faster 
than any polynomial in x. The idea is to multiply Eq. 
(EH) by and recall the relation il"la;"e'"-^ = d^e"^'''. 
Integration by parts and the inequality IjAip imply that 
the resulting integral converges and is bounded in x. This 
gives the decay. 

The polynomial symbols p{x,uj) = J2\^\=o 
with non-negative integer m correspond to differential 
operators of order m, p{x, dx) = X)™|=o 0'aix)dx- An ex- 
ample of a pseudodifferential operator that is not differ- 
ential is given by the symbol p{u;) — sin[ln(|a;|)], 
where fc is a real constant and is a cut function at 
\uj\ = 1/2. That is a smooth function that vanishes for 
\uj\ < 1/2 and is identically 1 for \uj\ > 1. The cut func- 
tion is needed to have a smooth function at w = 0. This 
symbol belongs to S'' . The function = x(ijj) ln(|tj|) 
is not a symbol, because \p{w)\ < co(w)'^, for every e > 0, 
but < ci((jj)~^, and the change in the decay is 

bigger than 1, which is the value of jaj in this case. An- 
other useful example to understand the symbol spaces 
is p{lli) — x(aj)|(jj|*'' ln(|ci;|), with k a real constant. This 
function is not a symbol for k natural or zero, for the 
same reason as in the previous example. However, it is 
a symbol for the remaining cases, belonging to S'^'^'^, for 
every e > 0. 

A very useful operator is A'* : S S given by 
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where s is any real constant. This is a pseudodifferential 
operator that is not differential. Its symbol is A* ~ (w) , 
which belongs to 5^, and then one says A* e V'*- It is 
usually denoted as A** = (1 — A)"/^. It can be extended 
to Sobolev spaces, that is, to an operator A^ : L^. 
This is done by noticing the bound ||A'*(m)|| = ||m||s for 
all M e 5, and recalling that 5 is a dense subset of L^. 
This operator gives a picture of what is meant by an s 
derivative, for s real. One can also rewrite the definition 
of H'^, saying that u £ H'^ ii and only if A'^{u) £ L^. 

Pseudodifferential operators can be extended to oper- 
ators acting on Sobolev spaces. Given p G tp"^, it defines 
an operator p{x,dx) ■ W . This is the reason 

to call m the order of the operator. The main idea of 
the proof is again to translate the basic estimate (|A1|I 
in the symbol into an L^-type estimate for the operator, 
and then use the density of S in L^. The translation 
is more complicated for a general pseudodifferential op- 
erator than for A* , because symbols can depend on x. 
Intermediate steps are needed, involving estimates on an 
integral representation of the symbol, called the kernel of 
the pseudodifferential operator. Pseudodifferential oper- 
ators can also be extended to act on distribution spaces 
5', the dual of Schwartz spaces S. 

An operator p : H^°° H^°° is called a smoothing 
operator if p{H~'°°) C C°° . That means p{u) is smooth 
regardless of u being smooth. One can check that a pseu- 
dodifferential operator whose symbol belongs to S~°° is a 
smoothing operator. For example, p{uj) — e"''^' S S~°°. 
However, not every smoothing operator is pseudodiffer- 
ential. For example, p{w) = p{^), with p G -ff" for some 
s and having compact support, is a smoothing opera- 
tor which is not pseudodifferential unless p is smooth. 
Friedrichs' moUifiers, Jc for e £ (0, 1], are a useful family 
of smoothing operators, which satisfy Jc(u) — > u in the 
sense, in the limit e ^ 0, for each u £ L^. 

Consider one more example, the operator X : S ^ S 
given by 

X{u) := [ e'"-^i|w|xHu(w)(?w, 

where xi'-^) is again a cut function at |a;| — 1/2. The 
symbol is A((jj) — i\uj\x{uj). The cut function x makes A 
smooth at Lu — 0. The operator without the cut function 
is ^ : iS ^ given by 

e{u) := / e''^-'^i\Lu\u{uj)dLU. 

Its symbol £{u!) — does not belong to any 5™ because 
it is not smooth at w = 0. Both operators A, i can be 
extended to maps L^. What is more important, 

their extensions are essentially the same, because they 
differ in a smoothing, although not pseudodifferential, 
operator. 

The asymptotic expansion of symbols is maybe the 
most useful notion related to pseudodifferential calcu- 
lus. Consider a decreasing sequence {'mj}°°^^, with 



limj^oo rrij = —oo. Let {p be a sequence of symbols 
p^{x,uj) £ S™^ . Assume that these symbols are asymp- 
totically homogeneous in co of degree rrij, that is, they 
satisfy p^{x,tuj) = t"''^ p^(x,uj) for |a;| > 1. Then, a sym- 
bol p £ 5*'"^ has the asymptotic expansion J2j Pj if ^^^d 
only if 

£-EEjj Vfe>l, (A3) 

and it is denoted by p ^ p_. . The first order term 
in the expansion, p^, is called the principal symbol. No- 
tice that nij are real constants, not necessarily integers. 
Every asymptotic expansion defines a symbol, that is, 
every function of the form '^^p . belongs to some sym- 
bol space S'"'^. However, not every symbol p £ 5™ 
has an asymptotic expansion. Consider the example 
p{u;) = x('^)|'^|"'^^^ lii(|cj|). The set of symbols that admit 
an asymptotic expansion of the form ljA3|l is called classi- 
cal, it is denoted by S"^, and the corresponding operators 
are said to belong to ^p'^. One then has C S"™. No- 
tice that if two symbols p and q have the same asymptotic 
expansion ^jP., then they differ in a pseudodifferential 
smoothing operator, because 

for all fc, and limfc^oo "t-j = — oo, so {p — q) £ S^°° . 
This is the precise meaning for the rough sentence, "what 
really matters is the asymptotic expansion." 

There is in the literature a more general concept of 
asymptotic expansion. It does not require that the p . to 
be asymptotically homogeneous. We do not consider this 
generalization in these notes. 

Most of the calculus of pseudodifferential operators 
consists of performing calculations with the highest or- 
der term in the asymptotic expansion and keeping careful 
track of the lower order terms. The symbol of a product 
of pseudodifferential operators is not the product of the 
individual symbols. Moreover, the former is difficult to 
compute. However, an asymptotic expansion can be ex- 
plicitly written for classical symbols, and one can check 
that the principal symbol of the product is equal to the 
product of the individual principal symbols. More pre- 
cisely, given p £ V'^j and q £ ip^v then the product is 
a well defined operator pq £ V'ci*^^ ^-^id the asymptotic 
expansion of its symbol is 

El^ -R^ [Wa;,w)] [a"g(a;,w)] . 

|a|>0 

Notice that the first term in the asymptotic expansion 
of a commutator [p, q\ = pq — qp, that is its principal 
symbol, is precisely 1/i times the Poisson bracket of their 
respective symbols, {p,g} = [dujopd^jq- d^jpd^jq). 
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Similarly, the symbol of the adjoint pseudodifferential 
operator is not the adjoint of the original symbol. How- 
ever, this is true for the principal symbols. The proof is 
based in an asymptotic expansion of the following equa- 
tion: 

J JR" 

There are three main generalizations of the theory 
of pseudodifferential operators present in the literature. 
First, the operators act on vector valued functions in- 
stead of on scalar functions. While this is straight- 
forward, the other generalizations are more involved. 
Second, the space of symbols is enlarged, first done 
in [23. It is denoted as S*^, and its elements satisfy 

\d^dZp{x,uj)\ < Cc,/3(w)"'-''l"l+*''^l, with Ca, ,0 a con- 
stant depending on the multi-indices a and (3. The extra 
indices have been tuned to balance two opposite tenden- 
cies; on the one hand, to preserve some properties of dif- 
ferential operators; on the other, hand to maximize the 
amount of new objects in the generalization. These sym- 
bol spaces contain functions like p(x,Lu) — {uj)°''^^\ which 
belongs to S^^^, where (5 > and rn = max^gRn a{x). 
Third, the domain of the functions p(x^Lo) is changed 
from R" X R" to X R", with f]~C R" any open 
set. A consequence in the change of the domain is that 
p : C^{Vi) C°°{^), so the domain and range of p are 
not the same, which makes it more difficult to define the 
product of pseudodifferential operators. These notes are 
intended to be applied to hyperbolic PDEs on M", which 
are going to be converted to pseudodifferential operators 
in S^, so there is no need to consider the last two gener- 
alizations. 



4. Further reading 



There is no main reference followed in these notes; how- 
ever, a good place to start is Notice that the notation 
is not precisely the one in that reference. The introduc- 
tion is good, and the definitions are clear. The proofs are 
d-ifficult to follow. More extended proofs can be found in 
j2i^ . together with some historical remarks. The whole 
subject is clearly written in ^|. It is not the most gen- 
eral theory of pseudodifferential operators, but it is close 
to these notes. A slightly different approach can be found 
in [sJl , and detailed calculations to find parametrices are 
given in '32'|. The introduction of l3t|| is very instructive. 
The first order reduction using A is due to Calderon in 
j2^ , and a clear summary of this reduction is given in 

The field of pseudodifferential operators grew out of a 
special class of integral operators called singular integral 
operators. Mikhlin in 1936 and Calderon and Zygmund 
in the begining of 1950s carried out the first investiga- 
tions. The field started to develop really fast after a 
suggestion by Peter Lax in 1963 'si'l, who introduced the 
Fourier transform to represent singular integral opera- 
tors in a different way. Finally, the work of Kohn and 
Nirenberg presented the pseudodifferential operators 
as they are known today, and they proved their main 
properties. They showed that singular integral operators 
are the particular case of pseudodifferential operators of 
order zero. Further enlargements of the theory were due 
to Lars Hormander [2?! l36l|. 
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(2§|. may clarify this. Consider the functions 

v{t,x) = sin(nt) sin(na;)/n*'^"'", 

w{t,x) = sinh(nt) sin(na;)/n*'^"'", 

with p > 1, n, constants, defined on t > 0, a; G [0, 1]. 
They are solutions of the Cauchy problem for wave equa- 
tion and the Laplace equation, respectively, with pre- 
cisely the same Cauchy data on t = 0, that is 

vtt-Vxx = 0, ^It^o = 0, Vt\t=o = sin{nx)/n'', 

wtt + W:,:^ = 0, w\t=<} = 0, 'Wt\t^o = s\a{nx) / nF . 

As n oo, the Cauchy data converges to zero in 
C^~^([0, 1]). In this limit, the solution of the wave equa- 
tion converges to zero, while the solution of the Laplace 
equation diverges. The concept of well posedness is in- 
troduced in order to capture this behavior of the wave 
equation's solution under high frequency perturbations 
on its Cauchy data. 



